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Abstract 

> ■ 

. The precise quark mass dependence of the one-loop effective action in an instanton background 

oo ! n 

has recently been computed 1]. The result interpolates smoothly between the previously known 

(N ■ 

extreme small and large mass limits. The computational method makes use of the fact that the 

in : 

single instanton background has radial symmetry, so that the computation can be reduced to 

^ '■ 

a sum over partial waves of logarithms of radial determinants, each of which can be computed 

P-i 

' numerically in an efficient manner. The bare sum over partial waves is divergent and must be 

43 : 



regulated and renormalized. In this paper we provide more details of this computation, including 



^ . both the renormalization procedure and the numerical approach. We conclude with comparisons 

H : 

of our precise numerical results with a simple interpolating function that connects the small and 
large mass limits, and with the leading order of the derivative expansion. 
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I. INTRODUCTION 



The computation of fermion determinants in nontrivial background fields is an important 
challenge for both continuum and lattice quantum field theory. Explicit analytic results are 
known only for very simple backgrounds, and are essentially all variations on the original 



work of Heisenberg and Euler 



For applications in quantum chromodynamics 



(QCD), an important class of background gauge fields are instanton fields, as these minimize 
the Euclidean gauge action within a given topological sector of the gauge field. Furthermore, 
instanton physics has many important phenomenological consequences j^ll, l| 9, 10]. Thus, 
we are led to consider the fermion determinant, and the associated one- loop effective action, 
for quarks of mass m in an instanton background. Here, no exact results are known for the 



mass 



lass depen 



full mass dependence, although several terms have been computed analytically in the small 

~ n - Q 



and large mass 



12L 1 1 31 ] limits. Recently, in |l( , the present authors presented 



a new computation which is numerical, but essentially exact, that evaluates the one-loop 
effective action in a single instanton background, for any value of the quark mass (and for 
arbitrary instanton size parameter). The result is fully consistent with the known small and 
large mass limits, and interpolates smoothly between these limits. This could be of interest 
for the extrapolation of lattice results [14], obtained at unphysically large quark masses, to 
lower physical masses, and for various instanton-based phenomenology. Our computational 
method is simple and efficient, and can be adapted to many other determinant computations 
in which the background is sufficiently symmetric so that the problem can be reduced to a 
product of one-dimensional radial determinants. While this is still a very restricted set of 
background field configurations, it contains many examples of interest, the single instanton 
being one of the most obvious. It is well known how to compute determinants of ordinary 
differential operators 00,0,0,0; but in higher- dimensional problems with partial 
differential operators, one must confront the renormalization problem since there are now an 
infinite number of 1-D determinants to deal with (even when the partial differential operator 
has a radial symmetry). 

In this paper we present more details of the results of . In Section |H] we define the 
renormalized effective action in the minimal subtraction scheme, as introduced by 't Hooft 
0, and _,e what is known about the s-naU and large mass lim it, In SeetionDJwe 
review how the single instanton background reduces the spectral problem to a set of radial 
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problems, and indicate how to regularize the effective action. This reduces the computation 
to two parts, one of which is analytic and the other is numerical. The analytic part con- 
cerns the renormalization of the effective action, and for this we use a WKB expansion as is 



developed in our earlier paper 20]. We stress that this renormalization computation, which 
constitutes Section HVI is analytic and exact, even though we use a WKB expansion, since 
we show that only the first two orders of the WKB expansion contribute. Section IVl presents 
details of the numerical part of the computation and shows how to combine the numerical 
part of the computation with the renormalization part to obtain the finite renormalized ef- 
fective action, which is plotted in Figure El In Section IVII we present a simple interpolating 
function that has been fit to our results, and we also compare our precise mass dependence 
with the mass dependence of the leading derivative expansion approximation, which was 
computed previously in j^J. The final Section contains some concluding comments. In Ap- 



pendix A we give some details (not fully given in [20]) which are needed in the approximate 
effective action calculation using the WKB phase-shift method. For this WKB analysis, the 
Schwinger proper-time framework |3j provides a natural way to implement the renormaliza- 
tion procedure consistently. Appendix B confirms that the same result is obtained in the 
regular and singular gauges for the instanton background. 



II. RENORMALIZED EFFECTIVE ACTION IN A SELF-DUAL BACKGROUND 

An instanton background field is self-dual, and self-dual gauge fields have the remarkable 
property that the Dirac and Klein-Gordon operators in such a background are isospectral; 
that is, they have identical spectra, apart from an extra degeneracy factor of four in the 
spinor case and zero modes present in the spinor case 0,1^ 22 |- Since the one-loop effective 
action is proportional to the logarithm of the determinant of the respective operator, this 
has the immediate consequence that it is sufficient to consider the scalar effective action 
to learn also about the corresponding fermionic effective action, for any mass value m. In 
particular, for a quark in a background instanton field, the renormalized one-loop effective 
action of a Dirac spinor field of mass m (and isospin |), r^ en (A;m), can be related to the 
corresponding scalar effective action, rf en (v4; m), for a complex scalar of mass m (and isospin 
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m 



Ti n (A; m) = -2 T> cn (A; m) - - In — 



(2.1) 



where p is the renormalization scale. The In term in (|2.1j) corresponds to the existence of a 
zero eigenvalue in the spectrum of the Dirac operator for a single instanton background. 

The one-loop effective action must be regularized. We choose Pauli-Villars regularization 
adapted to the Schwinger proper-time formalism, and later we relate this to dimensional 
regularization, as in the work of 't Hooft 61. The Pauli-Villars regularized one- loop scalar 

nn 

effective action is |12J, [13j 

V s (A- m) - In pet(-D 2 + m 2 ) Det(— <9 2 + A 2 ) j 

1 A (A, m) - m ^ Det( _ Q2 + m2) Det{ _ D2 + A2) j • {2.2) 

where D 2 = D^D^, with = <9 M - L4 M (x). In NT^ .A is a heavy regulator mass. We 
consider an SU(2) single instanton in the regular gauge 



A,{x) = A^xf- - 

_ a T a 



r 2 + p 2 ' 
2p 2 r] f , ua T a 



(r 2 + p 2 ) 2 ' 

where 77 Mi , are the standard 't Hooft symbols jfj 0] . 

The regularized effective action ()2.2j) has the proper-time representation 

r °° ds , 



(2.3) 



T b A (A;m) 



s 

00 ds 



e~ A2s ) I d A x tr(x|e- 



*(-d 2 ) _ e - s (~a 2 ) 



e~ m2s -e-^ s )F(s). 



(2.4) 



The renormalized effective action, in the minimal subtraction scheme, is defined as 



12 



rf en (A;m) 



lim 

A— >oo 



lim 

A— >oo 



T S A (A; rn)-±j±y 2 ln (^j J d 4 x tr(i^) 



ri(A;m)-^ln^ 



(2.5) 



where we have subtracted the charge renormalization counter-term, and p is the renormal- 
ization scale. By dimensional considerations, we may introduce the modified scalar effective 
action rf en (mp), which is a function of mp only, defined by 



rf en (A;m) 



r ren(™p) + g 



(2.6) 
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FIG. 1: Plot of the analytic small and large mass expansions for f^ n (m), from Equation (|2.7|) 
Note the gap in the region 0.5 < m < 1, in which the two expansions do not match up.. 

and concentrate on studying the mp dependence of T^ en (mp). Then there is no loss of 
generality in our setting the instanton scale p = 1 henceforth. 
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It is known from previous works that in the small mass (J 
limits, rf en (m) behaves as 

a ( |J + \ (In m + 7 — In 2) m 2 + . . . , m — > 



12j and large mass 



lnm 



1 



17 232 
+ 



7916 



(2.7) 



6 75m 2 735m 4 2835m 6 148225m 8 



+ 



m — > 00 



where 



a 



_A _ 2C'(-1) - - In 2 
72 s v ; 6 

0.145873... , 



(2.8) 



and 7 ~ 0.5772 ... is Euler's constant 24j. The leading behavior of the small mass limi; 
in (|2.7D was first computed by 't Hooft [6j, and the next corrections were computed in 11 1 
and [12]. This small mass expansion is based on the fact that the massless propagators in 
an instanton background are known in closed- form On the other hand, the large mass 
expansion in (|2.7|) can be computed in several ways. The 0{—^) and 0(=r) terms were 



computed in Q|, while the next two terms were computed in Q]. A very direct approach 
is to use the small- s behavior of F(s), the proper-time function appearing in (j2.4j) . as given 
by the Schwinger-DeWitt expansion. In our case this expansion reads 12 ] 



1 1 17 2 116 o 3958 , , . 

s^0+ : F(s) + —s + s 2 s 3 + s 4 + ---, 2.9 

K J 12 75 735 2835 44675 ' K J 

and using this series with ()2.4|) immediately leads to the large m expansion for T^ en (m) in 

(EIll- 
Equation ()2.7|) summarizes what is known analytically about the mass dependence of 
the renormalized one-loop effective action in an instanton background. This situation is 
represented in Figure^ which shows a distinct gap approximately in the region 0.5 < m < 1, 
where the small and large mass expansions do not match up. In this paper we present a 
technique which computes Tf en (m) numerically for any value of the mass m. Our results 
interpolate smoothly between the analytic small and large mass limits depicted in Figure 

III. RADIAL FORMULATION 

Our computational approach makes use of the fact that the single instanton background 
()2.3|) has radial symmetry jf]]. This has the important consequence that the computation of 
the regularized one-loop effective action ()2.4|) can be reduced to a sum over partial waves 
of logarithms of determinants of radial ordinary differential operators. Each such radial 
determinant can be computed by a simple numerical method, described in Section The 
physical challenge is to renormalize the (divergent) sum over partial waves. 

In the instanton background (J2.3)) . with scale p = 1, the Klein-Gordon operator — D 2 for 
isospin | particles can be cast in the radial form Q| 

where / = 0, ~, 1, §, • • • , and j = | l + \ |, and there is a degeneracy factor of (2/ + l)(2j + l) 
for each partial wave characterized by (Z, j)-values. [Note that 1(1 + 1) can be identified with 
the eigenvalue of L 2 = L a L a for L a = —^r]^ ua x^d u , and j(j + 1) with the eigenvalue of 
J 2 = (L a + T a )(L a + T a ) for T a = ^]. In the absence of the instanton background, the free 
operator is 

r Ffi in Aid 4. m 

(3.2) 



d 2 -> n^ e 



_ 3d_ 41(1 + 1) 



2 



6 



This radial decomposition means that we can express the Pauli-Villars regularized effec- 
tive action (|2.2|) also as 

*™) = ^.. (2 ' +1)(2 ' +2) p det l J +hdet Ug+ m2 

_ lndet f%ii> ±g) _ lndet } (3.3) 

Here we have combined the radial determinants for — I + |) and (/ + 1, j = (/ + |) — §), 
which have the common degeneracy factor (2Z + 1) (2Z + 2), so that the sum over / and j 
reduces to a single sum over I. In our actual analysis, as explained in detail below, we 
need to consider only a truncated sum over I with the expression ()3.3j) . and hence possible 
ambiguities as regards effecting the infinite sum over I become irrelevant. 

In Section we present a simple and efficient numerical technique for computing each of 
the radial determinants appearing in ()3.3|) . But to extract the renormalized effective action 
we need to be able to consider the A — > oo limit in conjunction with the infinite sum over 
This can be achieved as follows. Split the I sum in (13.3)1 into two parts as : 

L oo 

Ti{A-m) = £ Tl (l) (A;m)+ £ T s ul) (A;m) (3.4) 
i=o,i,... l=L +\ 

where L is a large but finite integer. In the first sum, which is finite, the cutoff A may be 

safely removed since for any given finite I ^ , Q] , 

lim det 

A^oo 

Mm det 1 J£T 7 a J = 1 ■ (3-5) 

Thus, the first sum in (13.4)) may be written without the regulator A as 

E rf(A-m) = (3.6) 
*=o,i,... 

l [ (U {1 l+ u +m 2 \ ( Hn+i n + m 

^ (21 + 1)(2I + 2) lndet jg±aL - + lndet J r + cc 2 '\ 

This sum can be computed numerically, and we find [see Section Ej that for any mass m 
it is quadratically divergent as L — > oo. This divergence is canceled by a divergence of the 
second sum in (|3.4|) in the large L limit, as we show in the next Section. 





IV. WKB ANALYSIS AND RENORMALIZATION 



In the second sum in (|3.4|) we cannot take the large L and large A limits blindly, as each 
leads to a divergence. However, we show in this Section that the A divergence is precisely 
of the counterterm form in the renormalized action in ([2.5)1 . and that the large L divergence 
is such that it precisely cancels the large L divergences from the large L limit of the sum in 
(13. 6 j) . Thus we obtain a finite renormalized effective action. 

The advantage of this technique is that the large A and large L divergences of the second 
sum in (13.4(1 can be computed analytically, using the WKB approximation for the corre- 
sponding determinants. The WKB approach to radial determinants was derived in (^(J up 
to third order in the WKB approximation, and the relevant results are reviewed below, and 
in the Appendix. It turns out that in the large A and large L limits we only need up to the 
second order in WKB. The large L limit of the second sum in (|3.4(1 can then be analyzed 
using the Euler-Maclaurin summation formula 2^|. 

It is convenient to express the second piece by the proper-time representation 

1 / 2 



E r l(o(^™) 



ds 



— e 



-A 2 ; 



Fr(s) 



(4.1) 



with 



/ 



Fr(s) 



dr 



\ 



E fM 

/,(s,r) = (2Z + l)(2i + 2)[/ (M+i) ( s ,r) + / 0+|)0 ( S ,r)]. 



/ 



(4.2) 
(4.3) 



Here the term J °° drfi(s,r), obviously related to Tf_^(A;m), can be found using the scat- 
tering phase shifts of the Schrodinger problem with the radial Hamiltonian in (|3.1(1 . Also, 
in considering the above infinite sum over partial wave terms, it is now crucial to have the 
contribution from partial-wave — I + |) and that from partial wave (I + |,Z) treated 
together as a package, as indicated in (|4.3(1 . For details on this, readers may consult the 
Appendix. We here only note that, for large I, the WKB approximation becomes exact and 
so in the large L limit we can use the WKB approximation (to an appropriate order) for 
Fl(s). According to the WKB expressions derived in the Appendix, it is found that, for 
each I, fi(s,r) has a local expansion in terms of the Langer-modified |27( potential, Vyj)(r), 



S 



and the corresponding Langer-modified free potential, V/(r) : 

4(/ + ±) 2 , 4(7 -00' + 1 + 1) 



r 2 + 1 



(r 2 + l) 2 



Specifically, the first order WKB result for fi(s,r) is obtained by using the form 



c(l) 



exp 



" s %i)( r ) 



2 A /7TS 



exp 



-sVKr) 



in the right hand side of (|4.3|) . and the second order WKB by using 

1 / s s 2 d 2 V(ij) x 



(4.4) 
(4.5) 

(4.6) 



2^ \4r 2 12 dr 2 
1_ / s s 2 d 2 V{ 



exp 



exp 



-sVAr) 



(4.7) 



20ri V 4r2 12 dr2 

The corresponding expression for the third order of WKB is also given in the Appendix, but 
this result is not needed for our present purposes. 

An important observation (which holds true to any order in the WKB expansion) is that 
the I dependence in f(ij)(s, r) has the form of a polynomial in / multiplied by an exponential 
in which I appears quadratically. Thus, if we use the the Euler-Maclaurin expansion [26 1 for 



00 r 00 1 1 

E fi = 2/ dlf(l)--f(L)--f'(L) + ... 



(4.8) 



all terms in this expansion, including the integral term, can be computed analytically. The 
integral term yields an error function 



dl exp 



-al 2 - 2bl 



^e^Erfcfe^ 
4a V 



(4.9) 



where Re(a) > 0. We here remark that it is important to perform the /-sum prior to 
considering the r-integration. For more details on the issue of the integration order, see the 
Appendix. 

To compute the sum in (j4.1j) we still need to perform the proper-time integral over s as 
well as the radial integral over r appearing in the WKB expression (|4.2jl . To achieve this, 
we adopt the following procedure. First, we trade the regulator mass A for a dimensional 
regularization parameter e, by demanding that 



fds[-l(e-^-e^)F Li s)}^fds[-l( 



S e )Fr(s) 



(4.10) 



Then, from the facts that Fl(s) ~ F(s) as s — > 0+, and F(s) = i + 0(s) for small s, we 
see that (jUTUj) requires 



6 \m 



-— + —h + 21nm) + 0(e). 
12e 12 u ; w 



(4.11) 



Thus, the correspondence between e and A is 



(4.12) 



7 + 21nA 

We now proceed to do the proper-time and radial integrals as follows. First, for L very 
large, it becomes convenient to rescale variables as s —>■ y/L 2 and r — > a;,/^. Second, we 
expand all terms (except the e~ m2s factor) in decreasing powers of large L. Then the y 
integral can be performed in closed form, yielding incomplete gamma functions [24]. These 
can be further expanded for large L, after which the x (that is, the radial) integral can be 
done. It is straightforward to perform these operations using Mathematica. To zeroth order 
in e, the results for the first and second order of WKB are given below. For the first order 
WKB term, the result is 



oo poo 

dr I ds 
o J o 



1 / 2 

— e ~ ms s £ 



J=L-\ 



1 „ AT hiL (\ 2 \ 119 ln2 V(l 

+ 2L 2 + AL - + m 2 + + 

24e 2 V6 / 72 12 24 

m 2 , , , 1 + 6 m 2 _ / 1 

+ _ (1 _ 2hl2 + lnm) ___ + 



where ^(5) — — 7 — 2 In 2. For the second order WKB term, the result is 



OO /"OO 

dr / 

Jo 



(e- m2s s e 



( 



(2), 



s,r 



V= 



1 InL 1 In 2 V'd! 



24e 



12 + 9 12 + 2 I 



12L + °VL 2 



(4.13) 



(4.14) 



The third order WKB term gives a contribution of at most O(-p-), and has no e pole. 
Similarly, it can be shown that all higher order WKB terms have no e pole, and vanish for 
large L. 

We can thus compute the large L limit by considering only the relevant parts of the first 
two WKB expressions in f!4.13|) and (|4.14|) . Inserting the identification between e and A in 
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(jHI21), we obtain 



m) 



-lnA + 2L 2 + 4L 
6 



1 ™U r 
- H In L 

.6 2 j 



+ 



127 1 m 2 

In 2 + — 

72 3 2 



m In 2 H In m 

2 



+ 



(4.15) 



We can now identify the physical role of the various terms in f)4.15j) . The first term is the 
expected logarithmic counterterm which is subtracted in ()2.5|) . and explains the origin of 
the |ln/i term in (j2.6|) . The next three terms give quadratic, linear and logarithmic diver- 
gences in L. We shall show in the next section that these divergences cancel corresponding 
divergences in the first sum in (|3.4j) . which were found in our numerical data. It is a highly 
nontrivial check on this WKB computation that these divergent terms have the correct co- 
efficients to cancel these divergences. Note that the InL coefficient, and the finite term, are 
mass dependent. 

Thus, the minimally subtracted renormalized effective action Tf en (m), defined in ()2.6|) . is 



lim J Y Tf(A;m) + 2L 2 + 4L- (- + — )\nL 



+ 



127 



1 , mr 
-ln2 + — 

3 2 



m 



m In 2 H In m 



(4.16) 



72 3 2 2 

where the first sum is to be computed numerically from the partial wave expansion in (j3.6|) . 



V. NUMERICAL CALCULATION 



In this section we describe the numerical technique for computing the radial determinants 
which enter the partial wave expansion in ()3.6j) . These one- dimensional determinants can be 



computed efficiently using the following result [15j, |16|, |17|, |18|, 1^|. Suppose M.\ and M. 2 are 
two second order ordinary differential operators on the interval r e [0, oo), with Dirichlet 
boundary conditions assumed. Then the ratio of the determinants of M.\ and M.2 is given 
by 



det 



lim 



\M 2 J r^\4> 2 (R), 
where ipi(r) (for i = 1,2) satisfies the initial value problem 

Mi^ir) = 0, with ^(0) = and ^'(0) = 1 



(5.1) 



(5.2) 
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Since an initial value problem is very simple to solve numerically, this theorem provides 
an efficient way to compute the determinant of an ordinary differential operator. Note in 
particular that no direct information about the spectrum (either bound or continuum states, 
or phase shifts) is required in order to compute the determinant. 

We can simplify the numerical computation further. Note that for the free massive Klein- 
Gordon partial- wave operator, T-&^ e + m 2 (with Ti { ^ e given in (|3.2|1 ). the solution to ()5.2|) is 
the modified Bessel function 



24] 



4r{r) = '^fl . (5 .3) 

This solution grows exponentially fast at large r, as do the numerical solutions to ()5.2|) for 
the operators 7~C(i t j) + m 2 , with Hn,j) specified in ()3.1|) . Thus, it is numerically better to 
consider the ODE satisfied by the ratio of the two functions 

This quantity has a finite value in the large r limit, which is just the ratio of the determinants 
as in (j5.1|) . The boundary conditions for the ratio function are 

n {lJ) (0) = l ; 7^)(0) = (5.5) 

A similar idea of considering the ratio function was used by Baacke and Lavrelashvili in 



their analysis of metastable vacuum decay 



21 



It is worthwhile making a brief side comment about the boundary conditions at r = 0. 
The two functions on the right hand side of ()5.4j) do not necessarily satisfy the boundary 
conditions at r = in (|5.2j) . For example, when I = 0, ^%ff )(r) does not vanish at r = 0. And 
in the massless case this issue is more serious. However, since only the ratio is important, 
one can introduce an ultraviolet regulator by imposing the boundary conditions in ()5.2|) at 
r = a, for a small but nonzero. Then the free solution satisfying the boundary conditions 
is a linear combination of the two modified Bessel functions l2l+1 ^ mr ^ and K2l+1 ^ mr ^ . It is 

r r 

straightforward to show that the differential equation governing the ratio function, and the 
asymptotic (r — > oo) value of the ratio function are independent of a as a — > 0. 

In fact, since we are ultimately interested in the logarithm of the determinant, it is more 
convenient (and more stable numerically) to consider the logarithm of the ratio, i.e., 

S { i tj) (r) =lnft ( y)(r) , (5.6) 
12 



which satisfies the differential equation 

^ + ( d Im)\(l + 2 Jjm^y^m = UM(r) , (5 . 7 ) 

dr 2 \ dr J \r i 2 m( mr )/ dr y ' 

with boundary conditions 

S(W)(r = 0) = , S' {lJ) (r = 0) = . (5.8) 

The 'potential' term U(j^{r) in (|5.7jl is given by 



4(j-/)(j + Z + 1) o 

f/ (^)( r ) = (r 2 + 1) 2 - ( 5 - 9 ) 



To illustrate the computational method, in Figure E]we plot S^ l+ i^{r) and ^(r) for 
various values of I, with mass value m = 1 (which is in the region in which neither the large 
nor small mass expansions is accurate). Note that the curves quickly reach an asymptotic 
large-r constant value, and also notice that the contributions from S^ ll+ i^(r = oo) and 
Sq+l = oo) almost cancel one another when summed. This behavior is generic for all 
values of mass m. 

To obtain very high precision for <SVj)(r = oo) in the numerical computation, it proves 
useful to make a further numerical modification. For large r, a good first approximation to 
S(i,j){ r ) 1S provided by neglecting the first two terms on the left-hand-side of the differential 
equation in (J5.7|) . Thus we define a new function Tn t j\(r) by 

= [ dr 'i^y) +T ^ ( 5 - 10 ) 

W l {r) = -+2m f +li ^ r \ . (5.11) 
r I 2 i + i{mr) 

This new function Tnj)(r) satisfies the modified equation 

d2T (i,j) | ( dT (u) \ 2 1 ( WM i o u dA r ) \ dT (i,j) _ ( u m( r ) \ 2 d (^Wr) / k i o\ 

with the boundary conditions: T(/j)(0) = T^^(O) = 0. Numerical values for the quantities 
defined in fl5.11j) provide greater accuracy at large r, and are also better for large /. In fact, 
we can iterate this type of transformation as many times as we wish. For our computation 
we achieved excellent numerical precision by iterating this transformation twice. 

The large r values of S^+i^r) and Sq + i ^(r) can be extracted with very good precision 
(we integrated out to r = 10 s ). Notice that the asymptotic values of Syj + i^(r) and S^ l+ i^(r) 
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FIG. 2: Plots of the r dependence of j + i^(r) and S^ + i ;)( r )> solutions of the nonlinear differential 
equation in l|5.7|) . for m = 1, and for / = 0, 10, 20, 30. The upper curves are for while the 

lower ones are for Sy + i ^(f). Note that the curves quickly reach an asymptotic large-r constant 
value, and also notice that the contributions from S^ l+ i^(r = oo) and S^ + i^(r = oo) almost 
cancel one another when summed. 

very nearly cancel one another, as illustrated in Figure El This behavior occurs for all m, 
and becomes more accurate as / increases. In fact, for a given mass, it is found that, as a 
function of /, S^^ l+ i^(r = oo) grows like In / while Sy + i^(r = oo) decreases like — In/. This 
divergence cancels in the sum, resulting in the behavior: 

P ( l ) = 5 (M+i)( r = °°) + *V±,/)( r = °°) 

~ O (\) , / -> oo (5.13) 



./, 

This behavior is illustrated in Figured which is a blow-up of the P(l) data in Figure El 

Recall from (jHIEj) that the first numerical sum in (|4.1fij) is in fact a sum over P(l), with 
degeneracy factors: 

E Tf{A;m)= £ (21 + 1)(2Z + 2)P{1) (5.14) 
;=o,i,... i=o,i,... 

Thus, we can rewrite our final expression ()4.16|) for the minimally subtracted renormalized 
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FIG. 3: Plot of Sq i + l-)(r = oo), S^ l+ i ^(r = oo), and their sum P(/), denned in ()5.13j) . for m = 1. 
Note that S^^i^r = oo) and <S*^ + i ^(r = oo) almost cancel, with their sum P(l) vanishing at 
large I. See also Figure 0J 
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FIG. 4: Plot of the I dependence of P(l), for m = 1. (This is a blow-up of the P(7) data from 
Figure EJ). P(l) behaves like O(j) for large I. Note that this implies that the sum over I in (|5.14|) 
has quadratic divergences in the large L limit. 



15 



effective action T^im) as 



(m) 



lim 

L— >oo 



{2l + l){2l + 2)P(l) + 2L 2 + AL 



1 m 2 \ , _ 
- H In L 

6 2 ; 



127 1 , m 2 

In 2 H 

72 3 2 



m 2 In 2 H In m 



(5.15) 



The first sum is over terms that are computed numerically, as described above. The rest 
represents renormalization terms which have been computed using minimal subtraction and 
WKB. 

Since the degeneracy factor (21 + 1)(2Z + 2) is quadratic, the large / behavior of P(l) 
indicated in (J5.13|) [and plotted in Figure 0] shows that in the large L limit, the sum (J5.14J1 
has potentially divergent terms going like L 2 , L and In L, as well as terms finite and vanishing 
for large L. Remarkably, we find that these divergent terms are exactly canceled by the 
divergent large L terms found in the previous section for the second sum in (|3.4|) . Thus, 
the renormalized effective action rf en (m) calculated using ()5.15j) is finite, converges for large 
L, and can be computed for any mass m. We found excellent convergence with L = 50 in 
our numerical data, combined with Richardson extrapolation 2^|. In Figure El we plot these 
results for rf en {m), and compare them with the analytic small and large mass expansions 
in (|2.7j) . The agreement is spectacular. Thus, our expression ()5.15|) provides a simple and 
numerically precise interpolation between the large mass and small mass regimes. 

As an interesting check, our formula (|5.15j) provides a very simple computation of 't 
Hooft's leading small mass result. In the m — > 0, ()5.7|) becomes 

2 



(1,3) 



+ 



dS (l,j) 



+ 



4Z + 3\ dS, 



dr 2 \ dr J \ r J dr 
One can find the solution of this equation in analytic form: 



U, 



(5.16) 



In 



2/ + 1 
2/ + 2 



+ In 



vi + 



21 + 1 vTTT 2 " 



In 



VT+ 



(5.17) 
(5.18) 



As r 



oo 



each quantity diverges but the sum, P(l), has a finite value 

'21 + 1 



P(l) = In 



2/ + 2 



(5.19) 
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FIG. 5: Plot of our numerical results for T^ cn (m) from ()5.15j) . compared with the analytic extreme 
small and large mass limits [dashed curves] from (|2,7|) . The dots denote numerical data points 
from (|5.15jl . and the solid line is a fit through these points. The agreement with the analytic small 
and large mass limits is very precise. 

Then it follows that 



f f en (m = 0) = Km | J2 (2/ + l)(2/ + 2)ln^|±l^+2L 2 + 4L-ilnL + ^-iln2 

17 1 1 
= _ii_ii n2 + i-2C , (-l) 
72 6 6 v ; 

= a f~) = 0.145873... (5.20) 
which agrees precisely with the leading term ()2.8|) in the small mass limit in (J2.7|) . 



VI. COMPARISON WITH OTHER RESULTS 

Since this is the first computation of the full mass dependence of the one-loop effective 
action in an instanton background, there is not much with which we can compare, except 
the small and large mass limits ([2.7)1 . which agree very well. There are, however, two other 
comparisons worth making. The first is with a modified Pade interpolating fit proposed in 
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12( , which is consistent with the two leading terms in each of the known analytic small and 



large mass limits given in (|2.7|l : 

~ c , , 1 \ mm + a — (3a + B)m 2 — 

riJm) nm + ^ —± ^ -S , 6.1 

rcnV 7 6 1 - 3m 2 + 20m 4 + 15m 6 ' V ; 

with a = a(l/2) ~ 0.145873 and (3 = |(ln2 - 7) ~ 0.05797. Based on the numerical 
data found in the present paper, we can fit the exact mass dependence in Figure El with an 
expanded form of this interpolating function. Let us assume the form 

~ s 1 l\nm + a-(3a + p)m 2 + A im 4 -A 2 m 6 

r„ n (m) ~ — lnm + - - . 6.2 

renV ; 6 1 - 3m 2 + Eim 4 + B 2 m 6 + B 3 m 8 K ' 

One may easily check that the leading two terms of the small mass expansion of this ex- 
pression ()6.2j) is the same as the small mass expansion in (j2.7|) . Then, comparing the four 
leading terms of the large mass expansion of (|6.2j) with the large mass expansion of (J2.7|) 
fixes the coefficients Bi,B 2 , B 3 to be: 

Z592955 255 „ 

B l = 25 A 2 + A x + 9a + 3(3 

V 21609 49 M 

B 2 = "75 ( ^-A 2 + A x ) , B 3 = 75A 2 . (6.3) 



,49 

There remain two free parameters, A\ and A 2 , unfixed in ()6.2j) . We can choose them so that 
the Pade approximant in ()6.2|) best fits the numerical data found in the previous section. 
This is a straightforward numerical exercise, and we find the best fit is given by 

A x = -13.4138, A 2 = 2.64587 (6.4) 

In Fig. |U1 we compare these approximations with the precise numerical data. Note that the 
fit based on ()6.2|) [solid line] is extremely precise, so we can use (16. 2|) as a simple analytic 
expression approximating the full mass dependence of the effective action, over the entire 
range of mass values. This is analogous to modified Pade fits used in chiral extrapolation of 
lattice data jl^ . and has also been explored for Heisenberg-Euler effective actions (2^. This 
form will also be useful if one wishes to capture the full p-dependence (for a given quark 
mass value) of the scalar effective action rf en (m) via (|2.6j) (and then also of the fermion 
effective action T^ en (m) via (j2.1|0 . 

Another comparison we can make is to the derivative expansion approximation. This 
approximation was already studied in [20], where it was noted that it was remarkably close 
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FIG. 6: Plot of "Pade" approximations of the effective action. The dotted line is the interpolation 
in (|6.1|) proposed in 2], while the solid line is the approximation in (|6,2j) which is an interpolating 
function fit to the exact mass dependence found in this paper. The solid dots are the exact 
numerical data. 



to the extreme small and large mass limits in (|2.7|) . Now that we have computed the 
full mass dependence of F"f en (m), it is worth revisiting this comparison. Recall that the 
philosophy of the derivative expansion is to compute the one-loop effective action for a 
covariantly constant background field, which can be done exactly, and then perturb around 
this constant background solution. The leading order derivative expansion for the effective 
action is obtained by first taking the (exact) expression for the effective Lagrangian in a 
covariantly constant background, substituting the space-time dependent background, and 
then integrating over space-time. For an instanton background, which is self-dual, we base 
our derivative expansion approximation on a covariantly constant and self-dual background 
1^ . This leads to the following simple integral representation for the leading derivative 
expansion approximation to the effective action 

/ , si 1 r°° dxx [ „ / 48a; 2 \ „ /-m 2 /4\/3x\ 

rf(A: m) = / <^ -84 + 14 In 1 + — - + 7V3 — arctan 



x 2 ( 7 11 48x 2 \l 1, 

+768— 2 Fr 1, -, -; r - - lnm 6.5 

m 4 \ 4 4 m 4 / 6 
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FIG. 7: Plot of f (m), comparing the leading derivative expansion approximation (solid line) with 
the precise numerical answers (dots). The dashed lines show the small and large mass limits from 



Figure [7| shows a comparison of this leading derivative expansion expression with the exact 
numerical data. In the range covered, the agreement is surprisingly good for such a crude 
approximation. 

VII. CONCLUDING REMARKS 

In this paper we have presented the details of a computation of the fermion determinant 
in an instanton background for all values of the quark mass. The agreement with the known 
analytic expressions in the small and large mass limits is excellent. As another application 
of our result, we can reinstate the dependence on p, the instanton scale parameter, simply 
by replacing m by mp. Then, given a quark field of fixed mass m, the fermion determinant 
as a function of instanton size p can be studied. For phenomenological applications, this can 
now be simply described with better than 1% numerical accuracy by using the interpolating 
function ()6.2j) for r^ n , together with our formulas (|2.6|) and (|2.1|) . The resulting interpolating 



(071) . 
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expression for the fermion determinant in an instanton background is 

e" r - n (7.1) 



rn 



— (pp) 3 exp < — -lnmp + 



1 \ \n(mp) + 2a- (6a + 2/3)(mp) 2 + 2A l [mpf - 2A 2 (mpf 



p K ^' r | 3 1 - 3(mp) 2 + B^mp) 4 + B 2 (mp) e + B 3 (mp) 8 )' 

This expression assumes minimal subtraction for the renormalized coupling entering the tree- 
level contribution. To obtain the corresponding expression for e~ r ™ n in other renormalization 
schemes, one needs to perform additional finite renormalizations, as discussed in Refs.(12, 
. Notice that such one-loop finite renormalization terms Arg nit0 can have dependence on 
the quark mass m (but not on p), and the lack of manifest decoupling for large quark mass 
in the expression ()7.1)1 is a renormalization artifact . We also remark that in instanton- 
based QCD phenomenology one may well choose the quark mass value m in ()7.1)1 to be 
different from the Lagrangian (or current) quark mass, taking instead some effective mass 
value [3^. Our formula ()7.1|) can be used for discussing instanton effects in gauge theories 



with compact extra dimensions as well 
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The computational method we described is versatile and can be adapted to a large class 
of previously insoluble computations of one-loop functional determinants in nontrivial back- 
grounds in various dimensions of space-time, as long as the spectral problem of the given 
system can be reduced to that of partial waves. One may especially consider using analogous 
methods for the computation of quantum corrections to the soliton energy in field theories. 
Several examples along this direction are currently under investigation. 
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APPENDIX A: ONE-LOOP EFFECTIVE ACTION BY THE WKB PHASE- 
SHIFT METHOD 



WKB theory is a powerful tool for obtaining a global approximation to the solution of a 



second-order ordinary differential equation 



Hence one expects that it can be utilized 



for the approximate calculation of a one-dimensional functional determinant |3J, |35[ . In the 
case of higher- dimensional functional determinants, which are usually needed in the one-loop 
effective action calculation of field theory, one can still try to use this WKB theory if the 
relevant partial differential operator becomes separable (as is often the case with rotationally- 
invariant background fields). In the latter case, however, no useful result can be derived 
from such analysis if one does not have an unambiguous renormalization procedure that 
goes with the WKB theory. In fact, the usual leading-order WKB theory is not sufficient for 
the determinant calculation if a consistent renormalization demands the contribution from 
higher-order WKB approximation to be included. We shall see below that this is the case. 

The renormalization problem mentioned above has been solved in our earlier paper j^j], 
by including needed higher-order WKB contributions within the Schwinger proper-time 
representation for the effective action. Also achieved there is a generalization of the 



Schwinger-DeWitt small proper-time expansion |32j,|33| to the appropriate expression for 
arbitrary proper-time value (in the case of rotationally-invariant background fields only) 
so that one can have an approximation to the full effective action. In this Appendix we 
shall briefly summarize this development and also provide further technical details on the 
formulas stated in j^J. 

In the proper-time representation ()2.4|) for the effective action, it is the function 

F(s) = J d 4 x tr(x| (e- s( - D2) - e~ s{ ~ d2) ) \x), (Al) 

which contains the important information. Given a rotationally invariant background field, 
we may utilize the phase shift analysis with scattering solutions of the 'Hamiltonian' 7i = 
— D 2 to rewrite the expression (jAl|) . To that end we will put the system in a large spherical 
box of radius R (with a Dirichlet or Neumann boundary condition at r = R), thus making 
the spectrum discrete (aj. In a single instanton background in particular, we may then 
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consider the quantum mechanical scattering solution for each partial wave, that is, 

(A2) 

The corresponding free Schrodinger equation yields 

= {"£ " f | + ^±2)Mt) = ^.M- (A3) 



We are interested in the solution of (jA2|) and (jA3|) that vanish as r 21 for r — > 0. Then, for 
large r, ipo(r) behaves as 

ip (r) ~ CV~ 3/2 cos[k (n)(r + a)], (A4) 

where 

ko (n + i)-k (n) = l + o(j^, (A5) 

(for nonnegative integer n), and a is a certain constant which is not important. On the other 
hand, we may write the large-r asymptotic behavior of the solution to (|A2|) in the form 

ip(r) ~ CV~ 3/2 cos[k(n)(r + a) + r){k{n))}, (A6) 

where r](k(n)) denotes the appropriate scattering phase shift [The related scattering matrix 
is give by S(k) = e 2M? ^]. Here, because of the boundary condition air — R, we may demand 
the discretized momentum k(n) to be related to ko(n) above according to 

k(n)(R + a) + r]{k{n)) = k {n)(R + a). (A7) 

From (|A7j) we conclude that 

M«) = fc(n) + ^+o(^)- (A8) 

If [k l 'i(n)} 2 and [&o(^)] 2 denote the energy eigenvalues associated with the Schrodinger 
equations (jA2|) and (jA3|) . respectively, the function F(s) give by (|Alj) can be represented as 

F ( s )= E E^ + ^^ + ^Ei 6 " 5 ^^ 12 - 6 "*^ 13 }' ( A9 ) 

i=o,i,— i n 
including the degeneracy factor (2/ + l)(2j + 1) for the partial wave. But, because of 
i8|). we find for large R 



e 1 v " — e 



- S [k^{n)? _ „-#)] 2 j = e - s [fc ! ^(n)] 2 f 2fc t -J(n)^-(fc(n)) j + j . ( A10 ) 
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Using (|A10|) in ()A9|) gives rise to 

F(s)~ E E( 2/ + l)(2j + l)E2(Afc)e- s[ " J(B)|2 2fc '' J( " )l? ' j(fc(B)) ^ (All) 

where AA; = -|, and then, replacing the sum J2 n by an integral for R — > oo, we obtain the 
following formula: 

2s 



- E E(^ + 1 )( 2 i + 1 ) / dke- k2s k Vld (k). (A12) 
* J=o,±... j 



With (|A12J) some caution must be exercised in dealing with the infinite partial wave 
sum. Actually, in the instanton background we are considering, the nature of the scattering 
problem as defined by (jA2j) and (jA3|) does not allow us to consider the /-sum and j-sum 
in (jA12|) in a completely independent manner. The point is that, as one looks at the given 
forms of H.{i,j) and 7i f ^ e , their small-r behaviors match for a given /-value; but it is the 
i-value that governs the large-, behavior of the potential entering ««,, while j does not 
appear in H { ^° at all. To obtain a convergent expression from (|A12J) . it is necessary |2(j] 
to consider the (/, j = I + |) and (/ + h,j = /) partial-wave contributions, both of which 
have the same degeneracy factor of (2/ + 1)(2/ + 2), together as one package. With this 
understanding, the expression (jA12|) can now be cast in the form 

2s r°° or i 

F(s) = - £ (21 + 1)(2Z + 2) dk e- k s k ki(fc) + Vi + ±J k )} ■ ( A13 ) 

If one has complete phase shifts for all partial waves at hand, one may use this formula 
()A13|) to calculate the function F(s) and then the one-loop effective action as well. But, in 
the massive case, the exact phase shifts cannot be obtained analytically. Therefore, in 20], 
we proposed to use the WKB expressions for the phase shifts, together with our formula 
(1A13I) . This method is elaborated below. 

n 

First, we need the results of Dunham [31] for higher-order WKB approximations to the 
scattering phase shifts. If the Schrodinger equation is written in the form 

I dx 



2 +Q(x)\*(x) = 0, (A14) 



the phase shift in the leading WKB approximation is given by 

1 



I sf0(x) dx - ('free') 



(A15) 
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where the integration path goes around the turning point r\ (i.e., the point where Q(x) 
vanishes) in the complex plane, and crosses the real axis at r = r$ (with ro < r) and 
r = r 2 (with r 2 taken to positive infinity). The choice of r has no effect on the value of 
the integral, and ('free') in (jA15|) represents the same integral but with Q(x) of the free 
Schrodinger equation. Dunham also derived the formulas in the second and third order 
WKB approximations: 

l Q"(x) 



,(2) 



V 



(3) 



48 Q(xf/ 2 
1 QW(x 



dx — ('free', 

7 [Q"{x)f 



dx 



f'free' 



(A16) 
(A17) 



768 Q{xfl 2 1536 Q{x) 7 / 2 
One cannot use Dunham's formula directly with the radial Schrodinger equation in (IA2I) 

_ n 

and (JA3j) . The latter should be transformed appropriately, following Langer 27]. Thus, 
writing r = e x and introducing the function 



\l/(x) = rip(r)\ 7 



e x ip(r = e x ), 



(A18) 



we recast (1A2J) as 
d 2 ^(x 



dx 2 



+ e 



2;r 



k - 



4(/ + i) 2 4(j _ /)(j + / + i) 



,2:r 



e 22: + 1 



+ 



e 2x + 1 



*(x) = 0. 



(A19) 



Since this is of the form (|A14|) . we can get the relevant scattering phase shifts simply by 
setting 

2 



Qd,j)(x) = e 2x Ik 2 



4(/ + |) 4 (j-/)(j + / + l) 



:,2x 



e 2x + 1 



(e 2x + I) 2 



(A20) 



and also, in connection with the free equation (|A3|) . 



Q z (x) = e 2l ' fc 2 



4(/ + | 



o2x. 



(A21) 



Then, using the original variable r = e x ' and integrating by parts, the desired WKB expres- 
sions for the phase shifts assume the form 
1 



%3 
(2) 



2 

1/1 1 

2 



fc 2 -V(y)(r) dr-('free'), 



1 



+ 



1 d'^Cr) 



1 



(A22) 
dr- ('free'), (A23) 
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128r 4 (p-y (y) ( r ))3/ 2 
7 ((PVnArW 1 



128r 2 



1536 V 
('free'), 



dr 2 



(P- W)) 7/2 



dr 2 (k 2 -V m (r)f/ 2 

1 g%j)(r) 1 
768 dr 4 (A; 2 - V m (r)fl 2 




where Vnj-\(r) is the so-called Langer modification 27] of the potential 



V, 



4(/ + |) 2 4 (j-/)(j+/ + l) 



,- + " r^+ i " ' ~ FTTF' (A25) 

and the term referred to ('free') denotes the integral expression appearing before but with 
V(ij)(r) replaced by 

4(/ + l) 2 



V,r) 



(A26) 



Note that Q { i,j){x) = r 2 [k 2 - V (lj) (r)} 



^,and Q°(x)=r 2 {k 2 - Vr(r)}|^. 
Using the results (fAl2jl - (fAl4j) in (fA"T3]) . we can simplify the expression by carrying out 
the /c-integration. First, consider the leading order WKB. The first contour integral in (|A22J) 
can be changed to the integral along the real axis over the interval (ri(k), oo), where r\{k) 
is a turning point, i.e., k 2 — Vn^\{r\{k)) = 0. On the other hand, fc-integral in (jA13|) 
runs from to oo. We may here change the order of integration, that is, perform the k- 
integral prior to considering the r-integration: in this case, the integration range for k would 
be over the interval (fci(r), oo), where k\{r) represents the value specified by the condition 
ki(r) 2 — V(ij)(r) = for a given value of r. Similar consideration may be given to the 
second contour integral in (|A23|) . the 'free' part. Then, observing that we obtain from the 
first contour integral 



(and similarly the form from the second contour integral), we are led to the following 

leading-order WKB expression for F(s): 

F {1) (s)= Y (2/ + l)(2/ + 2) [°° dr ^ 1 



k 2 - Vt 



-ldk-e-^k 
2 J IT 



k 2 - V, 



(A27) 



1=0,1 



/ ( ( j + i)(^r) + /J; |)0 ( S ./-; 



-^(i,i)( r ) 



sV,(r) 
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(A28) 
(A29) 



We use similar procedures to simplify the contributions coming from the second- and third- 
order WKB phase shifts in flA23|) and (jA24|) . For this, a particularly useful relation is 

1 * e -Wr ( -n + i )[ (K = M;V ., ) . (A30) 



7T 



fc 2 - V(r) 



n +2 7T 



As a result, we obtain the higher order WKB expressions for F{s), i.e., F®(s) and F^ 3 \s), 
which may be expressed again by the form (|A28J) but with 

'M)^' 1 ' ' 2^^s \4r 2 12 dr 2 J 2 v ^ v /i\4r 2 12 dr 2 J 1 J 

(3) e-'^)W [ 5s 2 _^_^yj(r) 7^ /^ (w) (r) ^ S 3 rf 4 %,)(r) 

T ^ S,r) l^^i |32r 4 48r 2 dr 2 

< f 5s 2 s 3 d 2 Vi(r) 7s 4 




2v / vt\/s 1 32r 4 48r 2 dr 2 1440 

In connection with using the formula (|A28j) . our discussion will not be complete without 
being clear about the order between executing the infinite series sum over / and performing 
the (improper) radial integral. This is a subtle point, and one possible way to settle the 
issue unambiguously would be to check explicitly which order gives rise to the known small- s 
behavior for the function F(s) correctly. As was asserted in [20], doing the /-sum before the 
r-integration yields the correct result. If instead one performs the r- integration first and then 
considers the /-sum, it gives a result differing from the correct small s expression of F(s) 
by -jk [Note that, although we used the WKB series for the calculation, the thus- found 
difference is an exact result since the order- dependent ambiguity is purely a high-energy 
phenomenon and the WKB series can be trusted in the high energy limit]. In view of this 
remark, the correct formula to be used in the WKB analysis of the effective action should 
read 

poo ^ r , 

F(s)= dr £ (2/ + l)(2/ + 2)[/ (M+|) ( S ,r) + / (/+ i /) ( S ,r)]. (A33) 

The result of using this formula for the fermion determinant in a single instanton background 
is presented in [20]. 
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APPENDIX B: THE INSTANTON DETERMINANT IN THE SINGULAR 
GAUGE 



In this Appendix we address the question of the gauge invariance of the determinant or 
the effective action. The proper time representation of the effective action in (J2.4j) is written 
in terms of covariant derivatives and is clearly gauge invariant. However, to formulate the 
partial wave expansion, as discussed in Section UTTl we chose a particular gauge (|2.3|) . Thus 
the partial wave expansion used in the main text does not possess manifest gauge invariance. 
In this Appendix we verify gauge independence by showing that we obtain precisely the same 
result in a different gauge. Since our computational method relies on the radial symmetry 
of the background field, we are restricted in which gauge we can choose. The choice in (|2.3j) 
is often called the "regular" gauge. But an instanton background also has radial symmetry 
in the so-called "singular" gauge: 



A?*(x) 



r 2 {r 2 + p 2 ) 



(Bl) 



in the singular gauge. In 



components with p or v equal to 4 



differs from rj^a in ()2.3j) only by the sign in the 



0, 0]. The gauge field in the singular gauge has 
singular behavior in the vicinity of r = 0. One may worry about the validity of the radial 
approach in the main text in the singular gauge because of this. However, it turns out that it 
is not so singular. Physically, the reason for this is the conformal invariance of the instanton 
background 0, |^ . The differential operator — D 2 in the instanton background (jBl|) in the 
singular gauge can be written as, setting p — 1, 

r _9 2 __3 9_ 4L 2 A(J 2 -L 2 -T 2 ) AT 2 



D = 

sing 



:b2) 



7" ^^7" 7" ^ f ^ ^ 7" ^ j 1 ^ ^ 2 ^ ^ 2 j ^ 2 

In the region of r ~ 0, the "potential" term in ()B2j) has 1/r 2 singular behavior. But the 
last two singular terms proportional to T 2 combine into T 2 /(r 2 + l) 2 , which is regular. We 
decompose l/(r 2 (r 2 + 1)) into 1/r 2 — l/(r 2 + 1); i.e., into singular and regular parts. The 
singular part, (J 2 — L 2 )/r 2 combines with the orbital part L 2 /r 2 . We can cast (jB2|) in the 
form 

4(L 2 



D 

sing 



d 2 



3d_ 4J 2 

r dr 



+ 



AT 2 



(B3) 



This should be compared with the corresponding operator in the regular gauge [see ()3.1|) ]: 

d 2 3d AL 2 A(J 2 - L 2 ) AT 2 
1 1_ 



D 



XV.il, 



r 2 + 1) (r 2 + If 
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Note that these two have the same form except that J 2 and L 2 are interchanged. Therefore, 
in the partial wave analysis the radial Hamiltonian (for isospin |) can be written as 

d 2 3 d | 4j(j + l) | 4(l-j)(j+Z + l) 3 



dr 2 r dr 



r 2 + 1 



1) 



(B5) 



There is a one-to-one correspondence between the eigenvalues for the sector (j, I) in the 
singular gauge and those for the sector (7, j) in the regular gauge, and they have the same 
multiplicities, (2j + 1)(2/ + 1). 

The regularized one-loop effective action in the singular gauge can be written as 



sing) TTlj 



(2/ + l)(2/ + 2) <Mndet 

I -i/sing 



sing 
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free 
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In det 
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Here we have combined the radial determinants for (I, j = I + |) and (/ + |, j 



(i+i) 



as in Section IIII1 and we have arranged the free determinants appropriately. But from the 
above arguments, we know that 



n 



njrcg 
rL (l',l)- 



(B7) 



Hence the singular gauge expression in ()B6j) is identical with the regular gauge expression 
in ()3.3j) . So the Pauli-Villars regularized one loop effective action has the same value in 
the singular gauge and in the regular gauge, and therefore the renormalized effective action 
()2.5|) . and the modified effective action (j2.6J) . has the same value in each gauge. 
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